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Evidence that the pseudogap (PG) in a 
near-optimally doped B^S^CaC^Og+a sample 
destroys the BCS logarithmic pairing instability 
[l] raises again the question of the role of the 
PG in the high-temperature superconducting 
cuprates [2]. The elimination of the BCS in¬ 
stability is consistent with the view that the 
PG competes with superconductivity. However, 
as noted in [1], the onset of superconductivity 
with a T c ~ 90 K suggests an alternative sce¬ 
nario in which the PG reflects the formation 
of short range pairing correlations. Here, we 
report results obtained from a dynamic cluster 
quantum Monte Carlo approximation (DCA) 
for a 2D Hubbard model and conclude that (1) 
the PG, like the superconductivity, arises due 
to short-range antiferromagnetic correlations 
and (2) contrary to the usual case in which the 
pairing instability arises from the Cooper insta¬ 
bility, here, the strength of the spin-fluctuations 
increases as the temperature decreases leading 
to the pairing instability. 

The superconducting transition temperature can be 
determined from the Bethe-Salpeter gap equation 

(k,U) n ,k',U>n')G(k',U n f)G(-k',-UJn’) (1) 

n'k' 

X (j>a(k ,CJ n ') = A a.(})(x{k' ) UJ n ). 

Here G{k,ui n ) is the dressed single particle Green’s func¬ 
tion, the irreducible particle-particle pairing vertex 
and k and u> n = (2 n + l)nT are the usual momentum 
and Matsubara frequencies, respectively. The tempera¬ 
ture at which the leading eigenvalue of Eq. (|T]) goes to 1 
gives T c and the corresponding eigenfunction <j> a (k,uj n ) 
determines the symmetry of the gap. In spin fluctua¬ 
tion theories the pairing vertex is approximated by an 
effective interaction 
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= ^ ^ X(V w m) (2) 

with the spin susceptibility and U a coupling 

strength. Various groups have used experimental data to 
model x(Qi u> m)i G{k,ui n ) and U in order to determine 
whether a spin-fluctuation pairing interaction is consis¬ 
tent with the observed T c values. 


Dahm et al. [3] used inelastic neutron scattering (INS) 
measurements for YBa 2 Cu 3 C> 6.6 to model the spin sus¬ 
ceptibility %(g, u! m ) and a one-loop self-energy approxi¬ 
mation to determine G. U was an adjustable parameter 
estimated from INS and ARPES data. Using the result¬ 
ing G and V e s in Eq. 0, they concluded that a spin- 
fluctuation interaction had sufficient strength to account 
for the observed T c . Nishiyama et al. @j used inelas¬ 
tic neutron scattering results for x(g,w) and solved the 
Eliashberg equations for the heavy fermion compounds 
CeCuSi 2 and Celrln 3 . For reasonable values of U, they 
found T c values which were again consistent with the 
notion that antiferromagnetic spin fluctuations were re¬ 
sponsible for pairing in these materials. In a recent paper 
Mishra, et al. [1. used angular resolved photoemission 
spectroscopy (ARPES) data for a slightly underdoped 
BSCCO (T c = 90K) sample to examine the effect of 
the pseudogap (PG) on the superconducting transition 
temperature and to determine whether a spin-fluctuation 
pairing mechanism could account for the observed T c . 
They found that the usual BCS logarithmic divergence 
associated with the propagators in Eq. 0 was destroyed 
by the pseudogap and the leading eigenvalue A d(T) re¬ 
mained small, and was essentially independent of temper¬ 
ature. This raises old questions regarding the interplay 
between the PG and superconductivity [2] which continue 
to be of interest [SHE]. Here, using the dynamic cluster 
approximation (DCA), we explore spin-fluctuation pair¬ 
ing in a Hubbard model which exhibits a PG. 

The two-dimensional Hubbard model we will consider 
has a near neighbor hopping t, a next near neighbor hop¬ 
ping t'/t = —0.15, an onsite Coulomb interaction U/t = 7 
and a filling (n) = 0.92. We will work in energy units 
where t = 1. The DCA calculations [3] were carried out 
on a 4 x 4 cluster and employed both continuous-time, 
auxiliary-field (CT-AUX) quantum Monte Carlo (QMC) 
[TUI and Hirsch-Fye (HF) QMC llj methods to solve 
the effective cluster problem m- In the DCA approxi¬ 
mation, where Ff r ? ( depends only on a finite set of cluster 
momenta I \, the fc-sum in Eq. ([lj gives m 

*%(K,w n ,K',w n ')^(K\u/ n )<l> a (K',w n ^ 


Here N c = 16 is the cluster size and the pairing kernel 
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G(fc, w n )G(—fc, — U! n ) has been coarse-grained (averaged) 
over the momenta k' of the DCA patches 

N ,_. 

X P 0 P (K^n) = ^Y, G ( K + k '^r l ) G (- K - k, ,-^n')- 
k' 

(4) 

For the parameters we have chosen, the uniform static 
susceptibility xid = 0, T) versus temperature, shown 
in Fig. E}b exhibits a peak at T* = 0.22 below which 




FIG. 1. Pairing in the presence of a pseudogap. (a) 

The uniform static spin susceptibility x(? = 0, T) versus tem¬ 
perature for (n) = 0.92 t' = —0.15 and U = 7 peaks at a 
temperature T* = 0.22 and decreases below this as the pseu¬ 
dogap opens, (b) The leading eigenvalue A a(T) of the particle- 
particle Bethe- Salpeter equation versus temperature (circles) 
from a DCA calculation of Tf r ^. The d- wave eigenvalue for 
the spin- fluctuation interaction with x(<7 5 w m ) the RPA spin 
susceptibility from Eq. [8] and U = 6.7 is shown as the solid 
squares. 

it decreases as T is reduced m- This behavior, seen 
in measurements of the magnetic susceptibility [ 15] and 
Knight shifts [16] of underdoped (hole) cuprates, reflects 
the opening of a pseudogap. ARPES experiments cans] 
find that this gap is anisotropic, opening in the antinodal 
regions of the Fermi surface. This behavior has also been 


seen in DCA calculations of the single-particle spectral 
weight [Hj 19]. In Fig. [TJd, the temperature dependence 
of the leading eigenvalue of the Bethe-Salpeter equation 
©> is shown as the circles. Its eigenfunction has d -wave 
symmetry and A d {T) approaches 1 at low temperatures. 
Thus this model system has a pseudogap that opens be¬ 
low T* and a d -wave eigenvalue that increases towards 1 
as T decreases. 

In addition to suppressing the q = 0 spin susceptibility, 
we find that the opening of the pseudogap destroys the 
low temperature BCS logarithmic divergence of the d- 
wave projection of the pairing kernel 

Pod{T) = ^ cj) d (K,uj n )xo P {K,u n )<j) d (K,u) n ) (5) 

C K,u> n 

Here, Xo P {K,uj n ) is defined in Eq. Q and (f> d (k,w n ) is 
the d -wave eigenfunction, which is approximated as 


4> d (k,uj n ) 


(cos k x — cos ky) \co n \ < J 
0 otherwise 


( 6 ) 


with J ~ 4 t 2 /U. A plot of Pod(T) versus T is shown in 
Fig- & and one can see that below T*, Po d [T) is sup¬ 
pressed as the pseudogap opens [2t)J [2T]. Here we have 
normalized Po d {T) to its value at a temperature T = 0.5f 
above T*. For comparison, the solid squares in Fig. © 
show Po d {T) for (n) = 0.85 which does not have a pseu¬ 
dogap and one sees the usual BCS logarithmic behavior 
(dashed curve). 

The absence of the BCS divergence in Po d {T) when 
there is a pseudogap is consistent with the finding of 
Mishra et al. [22] ■ However, as noted, they found that 
with this suppression, the spin-fluctuation pairing inter¬ 
action failed to give a superconducting transition. Based 
on this, they suggested that the pseudogap reflects the 
presence of short-range pairfield correlations which grow 
below T* and become coherent at T c . This behavior 
could be likened to the magnetic response of the large 
U half- filled Hubbard model. In this case, the forma¬ 
tion of local moments when the temperature drops below 
~ U/2 is seen in an increase in the expectation value of 
the square of the local moment ( S1) = ((| (n- 1 - — n^.)) ). 
In a similar way one can look for the onset of local pair 
formation as T decreases below the pseudogap temper¬ 
ature T*. Here with A\ +x/ = c] +xT c^ - c] +x4 cj t and 

A d = ( A l+x,e- A l+y,e + A l-x,£- A l-y,e)’ we have emu¬ 
lated (AjA d ) versus temperature. As shown in Fig. [ij), 
this correlation function does increase as the tempera¬ 
ture decreases. However, the four near neighbor pairfield 
correlations 

( A £+x,£ A £+x,£) = ~^{ n £ n £+x) — G(^£^£+x)i (?) 

contribute the dominant contribution to this increase as 
shown in Fig. ©>■ These results suggest that the pseudo¬ 
gap is more closely related to the formation of short range 
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FIG. 2. Destructure of the BCS logarithmic instability 
and nature of local pairing correlations, (a) The loga¬ 
rithmic BCS increase of the d- wave projection of the pairing 
kernel Pod(T) for (n) = 0.92 is suppressed by the opening of 
the pseudogap (circles). Here Pod(T) has been normalized to 1 
at a temperature T = 0.5 1 above T* . At temperatures below 
T* , where the pseudogap has opened, the BCS logarithmic 
divergence is suppressed. The solid squares show Pod(T) for 
a filling ( n) =0.85 where there is no pseudogap and one sees 
the usual logarithmic increase as the temperature decreases, 
(b) The temperature dependence of the local d- wave pairfield 
correlation function (AjJ'Ad) (circles). The observed increase 
in <A+A d > as T decreases below T* reflects the development 
of near neighbor AF correlations (squares). 


antiferromagnetic correlations than to local pair correla¬ 
tions in agreement with earlier ideas of Johnston [15] and 
more recent theoretical results 1301231. This identifica¬ 
tion of the PG with the development of short-range AF 
spin correlations is also consistent with the increase of 
the spin-susceptibility x(Q = (7 r,7r),w m = 0) as shown 
in Fig. [ 3 ] and as seen experimentally [23] , 

Returning to the question of whether the spin- 
fluctuation interaction, Eq. (J2|, can lead to supercon¬ 
ductivity when the logarithmic singularity of the BCS 
kernel is suppressed, we use DCA results for G{k,w n ) to 
construct 14ff(<Z,w m ). Here, following Mishra et al ., an 


FIG. 3. DCA spin susceptibility and RPA fit. The AF 

spin susceptibility xdca(Q = (7r, 7r ),ui m = 0) from the DCA 
calculation (circles) and the RPA fit, Eq. 086 , with U = 6.7 
(squares). The AF response continues to increase as T de¬ 
creases below T* leading to an increase of the spin- fluctu¬ 
ation interaction so that even though the BCS logarithmic 
increase of Po(T) is suppressed, the d-wave eigenvalue A<j(T) 
increases as seen in Fig. [2] 


RPA form for \ is used 


Xrpa(Q, w m ) 


X0_{Q,U m ) 

1 - UXo(Q,U m ) 


( 8 ) 


with 

P 

XoiQi^m) = w G(K + Q,u> n + u] m )G(K,Lo n ) , (9) 

c K 

where G(K,uj n ) = N c /NJ2k' G(K + k',u) n ) is the DCA 
coarse-grained Green’s function. The coupling U is esti¬ 
mated from the approximate fit of xrpa to xdca shown 
in Fig. [3] 

Then, replacing rf r ^ by V e g and using DCA Green’s 
functions, we solve the Bethe-Salpeter equation pi). Re¬ 
sults for A d(T) are shown (solid squares) in Fig.^2] We 
conclude that the increase in the strength of the pairing 
interaction V e g leads to an increasing A d{T) similar to 
that which is found using determined from the DCA 
calculation. Thus, in spite of the absence of the BCS log¬ 
arithmic increase in Pod(T), we find that the increase in 
the strength of the spin-fluctuations leads to an increase 
in A d(T) as the temperature is lowered. This differs from 
the results of reference [1J and we speculate that this dif¬ 
ference arises from a failure of their parametrization of 
G(k,tu n ) by ARPES data taken at 140 K as the temper¬ 
ature is lowered. 

To summarize, we have used DCA calculations for an 
under (hole) doped 2D Hubbard model, which exhibits a 
pseudogap, to see whether a spin-fluctuation interaction 
provides a reasonable approximation of the irreducible 
pairing interaction. In this calculation, the dynamic 
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mean-field cluster is such that charge density and strip¬ 
ing instabilities are suppressed, leaving antiferromagnetic 
and d -wave pairing as the dominant correlations. We 
find that while the pseudogap eliminates the usual BCS 
logarithmic divergence of the pairing kernel, a pairing 
instability arises from an increase in the strength of the 
spin-fluctuation interaction as the temperature decreases. 
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